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We study some dynamical properties of a classical time-dependent elliptical billiard. We consider 
periodically moving boundary and collisions between the particle and the boundary are assumed 
to be elastic. Our results confirm that although the static elliptical billiard is an integrable 
system, after to introduce time-dependent perturbation on the boundary the unlimited energy 
growth is observed. The behaviour of the average velocity is described using scaling arguments. 
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1. Introduction 

The idea of unlimited energy growth (also known as the Fermi acc eleration) du e to the repeated collisions 
of particles with a moving wall was introduced by Enrico Fermi [Fermi . 1949 ] in 1949 as an attempt to 
explain the acceleration of cosmic rays. He proposed a very simple model where charged particles could be 
accelerated by collisions with moving magnetic field structures. His original model was later modified and 
studied considering different approaches. Many of them take into account the inclusion of external fields, 
simplifications, damping coefficients, quantum and relativistic effects. 

One of the most important versions of this problem is the well known Fermi-Ulam Model (FUM). This 
model consists of a classical point particle of mass m, bouncing between two rigid walls. One of them is 
assumed to be fixed while the other one moves according to a periodic function. It is important to emphasize 
that this system has a very rich phase space structure in the sense that, depending on the initial conditions 
and control parameters, one can observe invari ant spanning curves, chaotic seas and Kolmogorov-Arnold- 
Moser (KAM) islands. Later on Pustylnikov Pustvlnikov . 11983. 1995|] rep l aced the fixed wall existing 



Holmesl . 



i\ rep i acea tne nxed. wail existmg 
1982 : lEversonl . Il986l : iLuna-Costa 



in the Fermi-Ulam model by a constant gra vitational field ^ 

1990l : iLeonel fc Livorati l2008l : IVincentl . I2OO0I ] the so-called bouncer model. Despite the similarity between 
the two models, there is a huge difference between them mainly regarding the average velocity of the 
particle for long time. In the Fermi-Ulam model it was rigorously proved by Pustylniko v that the ex i stence 
of invariant curves in the phase space always prevents the unlim ited energy growth IPustvlnikov . 1994 1 
although other work ers provided an evidence for such a conclusion Ulam . 196ll : Lichtenberg &: Lieberman , 
I972I : Douadvl . 19821 ]. On the other hand, for specific combinations of both control parameters and initial 



conditions the phenomenon of unlimited energy growth can be observed in the Pustylni kov bouncer model . 



This surprising result was later d iscussed and explained by Lichtenberg and Lieberman [Lichtenberg et al. 



1980l : Lichtenberg &: Lieberman , Il992] and can be easily understood by looking at the phase space. The 



1 



2 DIEGO F. M. OLIVEIRA,MARKO ROBNIK 



FUM has a set of invariant spanning curves limiting the size of the chaotic sea (as well as the particle's 
velocity), but such invariant tori are not observed in the bouncer model and the energy can typically grow 
unbounded. A natural extension of the one dimensional billiard models are the two-dimensional billiard 
systems. Basically they are classified (i) integrable, (ii) ergodic and (iii) mixed. In case (i) the phase space 
consists of invariant tori filling the entire phase space and typical examples are the circular and the elliptical 
billiard whose the integrability in the case of the circle comes from the angular momer itum conservation , 
and th e prod uct of the angular momenta with respect to the foci in case of ellipse [ Kamphorstl . 19991 : 



Koillei] . [l99i |. In case (ii) the time evolution of a single initial cori dition is enough to fill the phase sp ace 
and two examples are the Bunimovich stadium Bunimovichl . Il979 | and the Sinai billiard Sinail . Il970l]. In 
case (iii), the r e is a representative numbe r of b i lliards that present mixed phase space structure Saitol . 

One 



case (ni), tne r e is a representative numoe r oi Dnnaras tnat present mixed pnase space structure 
19821 : iRobnil Il983l : iLeonel k McChntockl . l2005l : iLooac et all . I2OOII . I2OO6I : lOliveira Leonel l201 0> 



important property in the mixed phase space is that chaotic seas are generall y surrounding Kolmogorov- 
Arno l d-Moser (KAM) is l ands which are confined by invariant spanning curves Lichtenberg &: Lieberman . 
19921 : IPustvlnikovl . Il98l Il995 ^. In particular such curves can cross the phase plane and partition it into 



several separated portions of the phase space. One of the main questions about two dimensional time- 
dependent systems is: Under what conditions the unlimited energy growth will be observ ed? In this sense, 



a conjecture was proposed by L oskutov-Ryabov-Ak i nshin (LRA) [Loskutov et all l200Cll | and later it was 



proved by Gelfreich and Turaev jCelfreich &: Turaev . 2008al b|. This conjecture, known as LRA-conjecture, 
states that the existence of a chaotic component in the phase space with static boundary is a sufficient 
condition to observe Fermi acceleration when a perturbation is introduced. Very recently Leonel and 
Bunimovich Leonel &: Bunimovich . 201Cll | extended the conjecture to the existence of a heteroclinic orbit 
in the phase space instead of the existence of a set with chaotic dy namics. Resu l ts tha t corroborate the 
validity of this conjecture include the time- dependent oval billi ard [Leonel et ai . 2009 1. stadium billiard 



Rvabov fc Loskutovl . [2010l | and Lorentz gas [Oliveira et Zl . l2010l |. Lenz et al. [Lenz et aLl . [2008l . l2009l . l2010l | 



observed that after a time-dependent perturbation is introduced the separatrix gives place to a chaotic layer 
and the particle's velocity can grow unbounded. Initia lly it was propose d that the acceleration exponent was 
controlled by the driven amplitude of the boundary [Lenz etal. . 2008]. Later, when extensive simulations 
were taken into account, the authors observed that when the number of collisions wit h the boundary i s 
large enough the acceleration expo nent is the same independ ent of the driving amplitude [ Lenz et 'al. , 2010l ] . 
Such result has been confirmed in fOliveira &: Robnild . l201ll | and it is studied in more details in the present 
paper. 

In this paper we revisit the problem of a classical particle confined in a time-dependent closed region 
of elliptical shape. Our main goal is to describe and understand the behaviour of the average velocity as a 
function of the number of collisions with the boundary using scaling arguments. 

The paper is organized as follows. In section [2] we describe the necessary details to define the four- 
dimensional mapping that describes the dynamics of the system, and our numerical results. Conclusions 
are drawn in section [3l 



2. The model and the map. 

In this section we discuss all the details required for the construction of the mapping. The two-dimensional 
driven elliptical billiard consists of a classical point particle of mass m confined into a closed region within 
which it is freely moving and is suffering elastic collisions with the time-dependent boundary. We stress 
that the particle is not affected by any external field and travels freely on a straight line until it reaches the 
boundary (see Fig. [1]). We describe the dynamics of the system in terms of a four-dimensional nonlinear 
mapping T(0„, a„, V^, tn) = (^'n+i, "n+i, Ki+i, ^n+i) that gives: the angular positioij^ of the particle On', 
the angle that the trajectory of the particle forms with the tangent line at the position of the collision 
a„; the absolute value of the particle velocity Vn = l^nl and the instant of the hit with the boundary 



^ 6 is the canonical parameter of the elhpse (called the eccentric anomaly in astronomy) and is not the angle tp between 
the position vector {x,y) and the x-axis (polar angle). The connection between 9 and the canonical parameter tp is tanip — 
(Bo Mo) tan 0. 
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Fig. 1. Illustration of five collision with the time-dependent boundary. The corresponding angles that describe the dynamics 
are also illustrated. 



Leonel et al . 2009l |. The index n denotes the n^^ coUision with the moving boundary. The Cartesian 



components of the boundary at the angular position {6n,tn) are 



X{en, tn) = [Ao + C sin(t„)] cos(^„) , (1) 
) = [Bo + C sin(t„)] sin(6'„) , (2) 

where and Bq are constants, thus, at any time t„ we have elhptical shape. The control (oscillation) 
parameter < C < min^Ao, Bq) controls the amplitude of oscillation and G [0, 27r) is a counterclockwise 
canonical angle measured with respect to the positive horizontal axis. The angle between the tangent of 
the boundary at the position {X{9n),Y{9n)) measured with respect to the horizontal line is 



arctan 



Y'i 



x'{e„) 

where the expressions for both X'(9n) and Y'(9n) are written as 



(3) 



X\en) = -[Ao + Csm{tn)] sin(e„), (4) 
ViOn) = +[Bo + Csm{tn)] cos(^„). (5) 

Since both 0„ and an are already known, the angle between the horizontal axis and the particle's 
trajectory is {4>n + On) and the vector velocity of the particle is given by 

l^n = I Vn\[cos{<pn + an)i + sin(</)„ + an)j] , (6) 

where i and j represent the unit vectors with respect to the X and Y axis, respectively. This allows us to 
obtain the position of the particle as function of time as follows 

n 

I COs{<pn + an){t - tn) , (7) 
Yp(t) =Y{9n,tn) + \fn\sm{<Pn + an){t-tn) , (8) 

where the index p denotes that such coordinates correspond to the particle. In order to find the angular 
position 9n+i at the collision n + 1, and the time tn+i of the next collision, we need to solve numerically 
the expression Rp{9n+i,tn+i) = Rb{9n+i,tn+i) where Rp{9n,tn) is the distance of the particle measured 



with respect to the origin of the coordinate system, i.e., Rp{9n,tn) = \ Xp{t) + Yp{t), and by evaluating 
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Fig. 2. (Color online) Behaviour of Vvs.n for different initial velocities. The control parameters used were Aq — 2 , Bq = 1. 
the expression 

V[Xp{tn+l) - X{9n,tnW + [^^(Wl) " >^(^n,tn)]2 



n+1 



tn + 



This is equivalent to solving the two equations Xp{tn+i) = X(6n+i,tn+i) and Yp{tn+i) = Y{9n+i,tn+i)- In 
the comoving frame of the boundary, where the velocity of the particle for the elastic collision is denoted 
by 1^', the following conditions must be satisfied 



n--'-^ n+1 



At the new angular position, 9n+i, the unitary tangent and normal vectors are 



— > 



Tn+i = cos(0„+i)i + sin(0„+i)j , 
^n+i = -sin(0„4.i)i + cos(0„+i)j , 
and after some algebra we can easily find 

"l^n+l-^n+l = I I [C0S((?:>„,+1 -(/>„- a„)], 
"l^n+l-^n+l = I I [sin((/)„+i - - On)] + 2}^ h.J^ n+1, 

where is the velocity of the boundary which is written as 

'Pfj = Ccos(tn+i)[cos(6'„+i)i + sin(6'„+i)j] . 
Thus, the absolute value of the velocity just after the collision (n + 1) is given by 

Vn+l =1 \^n+l 1= \/[fn+l.'^n+l? + [K+l.i^n+l]^ , 

and the angle a„+i is 

n+1 

tana„+i = ^ — . 

V n+l-T n+1 



(10) 

(11) 

(12) 
(13) 

(14) 
(15) 

(16) 
(17) 

(18) 
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Fig. 3. Behaviour of V x n for different initial velocities. The control parameters used were Aq = 2, Bq = 1, C = 0.2. 



2.1. Numerical results 

Our numerical results for the time-dependent driven elliptical billiard discuss basically the behaviour of 
the average velocity of the particle as a function of the number of bounce^ n. This is done because in a 
4D phase space it is difficult or even meaningless to consider individual pointwise initial conditions. We 
obtain the average velocity in two steps. Firstly, we evaluate the average velocity over the orbit for a single 
initial condition 

n 

^^ = ;r-TE^M' (19) 

j=0 

where the index i corresponds to a member of an ensemble of initial conditions. Second, we take the average 
over the ensemble of initial conditions, so that the average velocity is defined as 

_ 1 *^ 

^=mE^- (20) 

i=l 

where M denotes the number of different initial conditions. We have considered M = 200 in our simulations 
and fixed the control parameters as Aq = 2, Bq = 1. We will concentrate on the infiuence of the initial 
velocity Vq on the behaviour of the average velocity. Our main goal is to describe such a behaviour using 
scaling arguments. Figure [2] shows the behavior of the average velocity as a function of the number of 
collisions for different initial conditions and different values of the control parameter C. Figure [3] shows the 
behavior of the average velocity V as a function of the number of collisions n for different initial velocities 
and fixed C = 0.2. We have chosen 11 different values for initial velocity Vq while a random choice for the 
other variables was made as f G [0, 27r], 6 £ [0, 2tt] and a € [0, vr]. As one can see, all curves of the V behave 
quite similarly in the sense that: (a) for short n, up to n ~ n^; the average velocity remains constant equal 
to Vip, where ip means initial plateau, and (b) after a crossover for n » n^, all the curves start to grow 
with the same exponent. For such a behaviour we propose the following hypotheses, which turn out to 
describe the empirical facts: 



■^It should be noted that as long as the dynamics is chaotic there is not fundamental difference in choosing the discrete time n 
or the continuous time t. Namely, if oc , then V oc t^/^^~^^ that is to say the acceleration exponent changes /3 — >■ /3/(l — /?). 
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Fig. 4. (a) Plot of Vip x Vg. (b) Behaviour of Ux as function of Vq. 



(21) 



(22) 



(1) For short n, say n <C rix, V behaves according to 

Vip cx Vq' , 

and if the initial plateau is weU defined we of course expect 7 = 1. 

(2) For Ux, the average velocity is given by 

V (xn^ , 

where the exponents 7 and /3 are critical exponents, namely the initial and acceleration exponents, 
respectively. 

(3) The crossover iteration number n^, that marks the change from constant velocity at the initial plateau 
to the growth is written as 

Tlx oc , (23) 
where z is the third exponent, called the crossover exponent. 



With these three assumptions, from the method of [Leonel. McClintock &: da Silval . l2004l | . we suppose that 
the average velocity is described in terms of a scaling function of the type 

y(yo,n) = Ay(A"Vb,AM , (24) 
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Fig. 5. (a) Behaviour of the average velocity for different values of Vq, a subset of the curves of figured (b) their collapse 
onto a single and universal plot. 

where A is the scahng factor, a and b are scaling exponents. If we chose A^Vq = 1, then X = Vq and Eq. 
is given by: 

V{Vo,n) = V-'^''V,{V-'K) , (25) 

where = V^^' 

'n) is assumed to be constant for n <^ Hx- Comparing Eq. (j25p and Eq. 
OTI) . we obtain 7 = —1/a. On the other hand, if we chose X^n = 1, which means A = n^^/^ and Eq. (p^ 
is rewritten as 

V{Vo,n) = n-^/'V2{n-^/''Vo) , (26) 

where the function V2 is defined as V2{n~'^^^Vo) = y(n~"/^Vo, 1). It is assumed to be constant for n ^> n^. 
Comparing Eq. (j26p and Eq. (I22p we find /3 = —1/6. Considering the two different expressions for the 
scaling factor A, the condition of the crossover exponent is A = n^; = fix =Vq ^^"^ = ^ and we find 



7 



(27) 



Note that the scaling exponents are determined if the critical exponents 7 and f3 are numerically 
obtained. The exponent /3 is obtained from a power law fitting for the average velocity when n ^ Ux- 
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Thus, an average of these values gives /3 = 0.50(1) (see Fig. [3]), and j3 = 0.49970(2) for others values of 
C (Fig. [2]). Figure [4] shows the behaviour of (a), VipVS.VQ and (b), UxVS.Vq. From a power law fitting we 
obtain z = 2.0123(1) = 2 and 7 = 0.99991(8). Considering the previous values of both /3 and 7 and using 
Eg. 1271 we find that z = 1.99(4). This result indeed agrees with our numerical data. A confirmation of the 
initial hypotheses comes from a collapse of different curves of Vvs.n onto a single and universal plot, as 
demonstrated in Fig. [5l Additionally, considering the fact that the critical exponents are /3 = 0.5, 7 = 1 
and z = 2 we can conclude that the conservative time dependent elliptical billiard belong s to the same 
class of universality like the conservative time-dependent Lorentz gas [Qliveira et al . 2010l ] for the range 
of the control parameters considered. 



3. Final Remarks 

As a conclusion of the present paper, we have studied some dynamical properties of a time-dependent 
elliptical billiard considering elastic collisions with the boundary. We have shown that the average velocity 
for short time remains constant but after a crossover it starts to grow with exponent 0.5. We have shown 
that such a behaviour can be described using scaling properties and we have found an analytical relation 
between the critical exponents /3, 7, and z (acceleration, initial and crossover exponents, respectively). 
Our scaling hypotheses are confirmed by a good collapse of all the curves of the average velocity onto a 
single universal plot, therefore confirming that the model is scaling invariant and also that it belongs to 
the same class of universality of the time-dependent Lorentz gas for the range of the control parameters 
consider ed. It should be noted that the physical origin and explanation of the scaling laws discovered 
by [ Leonel. McClinto ck fc da Silval '2004*1 is still largely unknown, although there is some progress on the 



theoretical side Batistic fc Robnik . .2011]. 
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